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ABSTRACT 

We study resolution effects in numerical simulations of gas-rich and gas-poor major merg- 
ers, and show that the formation of slowly-rotating elliptical galaxies often requires a resolu- 
tion that is beyond the present-day standards to be properly modelled. Our sample of equal- 
mass merger models encompasses various masses and spatial resolutions, ranging from about 
200 pc and 10 5 particles per component (stars, gas and dark matter), i.e. a gas mass resolution 
of ~ 10 5 Mq, typical of some recently published major merger simulations, to up to 32 pc 
and ~ 10 3 Mq in simulations using 2.4 x 10 7 collisionless particles and 1.2 x 10 7 gas par- 
ticles, among the highest resolutions reached so far for gas-rich major merger of massive disc 
galaxies. We find that the formation of fast-rotating early-type galaxies, that are flattened by a 
significant residual rotation, is overall correctly reproduced at all such resolutions. However, 
the formation of slow-rotating early-type galaxies, which have a low residual angular mo- 
mentum and are supported mostly by anisotropic velocity dispersions, is strongly resolution- 
dependent. The evacuation of angular momentum from the main stellar body is largely missed 
at standard resolution, and systems that should be slow rotators are then found to be fast rota- 
tors. The effect is most important for gas-rich mergers, but is also witnessed in mergers with 
an absent or modest gas component (0-10% in mass). The effect is robust with respect to our 
initial conditions and interaction orbits, and originates in the physical treatment of the relax- 
ation process during the coalescence of the galaxies. Our findings show that a high-enough 
resolution is required to accurately model the global properties of merger remnants and the 
evolution of their angular momentum. The role of gas-rich mergers of spiral galaxies in the 
formation of slow-rotating ellipticals may therefore have been underestimated. Moreover, the 
effect of gas in a galaxy merger is not limited to helping the survival/rebuilding of rotating 
disc components: at high resolution, gas actively participates in the relaxation process and the 
formation of slowly-rotating stellar systems. 
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1 INTRODUCTION 

Numerical simulations have been intensively used for more than 
two decades to study the properties of the remnants of galaxy 
mergers and the role of hierarchical merging in the formation 
of elliptical-l i ke early-type galax ies dHern quist & B arnes! Il99ll ; 
IBarnesI 1 19921 : IMihos et all 1 19951) . With the increasing resolu- 
tion and large statistical samp les (e.g.. INaab & Burkertl |2003| ; 
iBournaud. Jog & Combesl2005l : ldi Matteo et al.l200ll2008l) . mod- 
ern work tends to quantify in details the properties of major and 
minor merger remnants, and accurate comparisons with observed 
properties of ear ly-type galaxies can now be envisioned (e.g., 
IBurkert et alj2008[) . 

A general concern, though, is that the impact of the spatial 
and the mass resolutions on the detailed properties of the systems 
under scrutiny remains largely overlooked, and whether or not sim- 
ulations of mergers have converged with today's typical resolution 
remains unexplored. Obviously, increasing resolution enables sim- 
ulations to directly resolve cold gas clouds and clustered star for- 
matio n (e.g. lBour naud. Due & Emsellem 2008; Kim. Wise & Abell 
2009), but whether these additional small-scale ingredients can sig- 
nificantly impact the global, large-scale morphology and kinemat- 
ics of merger remnants has not been studied in detail. In cosmolog- 
ical simulations, an increase in resolution (i.e., an increase in the 
number of particles and/or decrease of the softening length) can 
affect the baryonic density and circular ve l ocity profiles of indi- 
vidual galaxies in a halo Waab et al.ll2007t) . iNavarro et ail d2008h 
also studied numerical convergence via a suite of ACDM simu- 
lations and confirmed that the halo mass distributions were better 
described by Einasto profiles that are not, stricly speaking, univer- 
sal. 

While many resolution studies have been made in cosmologi- 
cal simulations, f ew have focused on g alaxy mer ger simulations. 
ICoxetalJ j2006t) . iHopkins et ail J2008I) and e.g.. Idi Matteo et all 
d2008h included some checks of the effect of resolution on the 
star-formation activity of on-going mergers. But a resolution study 
aimed at examining the detailed morphology and kinematics of re- 
laxed merger remnants (i.e., galaxies which tend to be roughly S0- 
or elliptical-like) has not yet been conducted. 

Models of galaxy mergers have reached particu l arly h igh res- 
olution with the work of IWetzstein. Naab & Burk ert ( 2007) (70 pc 
softening with 4 x 10 6 particles in total - but o nly 45,000 for the 
gas component). iLi. Mac Low & Klessenl fc004l) (1 to 100 pc and 
5 x 10 5 gas particles per galaxv). lNaab et alj J2007l) (8 x 10 6 parti- 
cles with a 125 pc resolution in a cosmological re-simulation of an 
individual galaxy halo). The highest r esolution for gas-rich merg- 
ers ha ve been achieved recently by IBournaud. Due & Emselleml 
( 2008) for mergers of bright spiral galaxies, with a total of 3.6 x 10 7 
particles includin g more than 10 7 gas p article, and a 32 pc soften- 
ing size, and iKim. Wise & Abell 1 20091) with a spatial resolution 
of 3.8 pc and a mass resolution of 2 x 10 3 Mq (for dwarf or 
low-mass spirals, though). But such high-resolution studies have 
focused on small-scale gas physics and structure formation, with- 
out studying the impact of high resolution on the global properties 
of the elliptical-like galaxies formed in major mergers. 

Large samples of simulations of idealized galaxy mergers re- 
main typically limited to softening lengths of ab out 100-300 pc, 
and ~ 10 particles per galax y (see s amples in 
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di Mat teo et al.l2007T) . Whether or not the relatively limited numer- 



ical resolution used in such studies affects the global properties of 



merger remnants is still a largely open question: for instance, the 
detailed comparison of major m erger remnants with the observed 
anisotropy-flattening relation by Burk ert et alj d2008h relies on sim- 
ulations with 2 x 10 4 gas particles per galaxy, a gas particle mass 
~ 3 x 10 J Mq, and a spatial resolution (softening) of about 200 pc. 

Within the context of the ATLAS 30 project 
(http://purl.org/atlas3dl, an extensive set of numerical simu- 
lations is being conducted to support the multi-wavelength survey 
of a complete sample of early-type galaxies within the local 
(40Mpc) volume, in terms of various formation mechanisms 
of early-type galaxies: binary mergers, multiple mergers, disc 
instabilities, etc. An ambitious series of simulations of mergers 
are being specifically performed and analysed for this purpose 
(Bois et al. in preparation). To properly interpret the results from 
these simulation efforts, as well as to understand the robustness 
of existing and past studies of galaxy mergers, we first probe the 
effect of spatial and mass resolutions on the global structure of 
binary disc merger remnants. 

In this paper, we study the effect of numerical resolution on the 
global morphology and the kinematics of the simulated remnants 
of binary, equal-mass major mergers. We wish to examine resolu- 
tions ranging from the typical resolutions used in recent, large sim- 
ulations samples, to some of the highest merger simulations ever 
performed. We study both dry (collisionless) and wet (gas-rich) 
mergers of disc galaxies. The modeled interaction orbits lead to the 
formation of both fast rotators, i.e., early-type galaxies flattened 
by significant rotational support, and slow rotators, i.e., early-type 
galaxies with low residual rotation, supported (and flattened) by 
(anisotro pic) velocity dispersion s, following the classification de- 
tailed in Ems ellem et alj d2007l . see also Section [2j. We find that 
the formation of fast rotators is overall correctly reproduced with 
numerical simulations at modest resolutions. In contrast, the forma- 
tion of slow-rotating systems is correctly reproduced only at high 
resolution (Section|3j, above the resolution of most of the recently 
published merger simulations. The influence of gas on the struc- 
ture of merger remnants, compared to dry mergers, also differs at 
high resolution, and is not limited to easing the survival and/or re- 
building of rotating disc components. In Section|4] we further ex- 
amine the origin of this observed resolution effect in the formation 
of slow-rotating systems. We show that it is not an artefact from 
different initial conditions or interaction orbits, but that the phys- 
ical treatment of the merging process is actually biased when the 
resolution is too low. The effect of the resolution has been tested 
on other simulations producing slow rotators and we find that it is 
a systematic one (Section |5). We summarize our results, discuss 
the required resolution for accurate studies and the general impli- 
cations for the formation of elliptical galaxies in Section|6] 



2 SIMULATIONS AND ANALYSIS 
2.1 Method 

2.1.1 Code 

We use the particle-mes h code described in 
IBournaud. Due & Emselleml J2008) and references therein. 

This code uses a Cartesian grid on which the particles are 
meshed with a "Cloud-In-Cell" interpolation. The gravitational po- 
tential is computed with an FFT-based Poisson solver and particle 
motions are integrated with a leap-frog algorithm, and a time-step 
of 0.5 Myr. 

Interstellar gas dynamics is modelled with the sticky-particle 
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scheme with elasticity parameters fit = fir — 0.6. This scheme 
neglects the temperature and thermal pressure of the gas, assuming 
it is dominated by its turbulent pressure, which is the case for the 
star-forming interstellar medium at the scales that are studied here 
dElmegreen & Sc alo 2004; Burkert 2006). The velocity dispersion 
of the particles model the turbulence, and their mutual collisions 
are inelastic to ensure that the turbule nce dissipates over about a 
vertical crossing time dMacLowlll999ri . 

The star formation rate is computed using a Schmidt- 
Kennicutt law: it is then proportional to the gas density in each 
cell to the exponent 1.5. Gas particles are converted to star par- 
ticles with a corresponding rate in each cell. Energy feedback 
from supernovae is accou nted for with the scheme proposed by 
Miho s & Hernquistl dl994l) . Each stellar particle formed has a num- 
ber of supernovae computed from the fraction of stars above 8 Mq 
in a Miller-Scalo IMF. A fraction e of the 10 51 erg energy of each 
supernova is released in the form of radial velocity kicks applied 
to gas particles within the closest cells. We use e = 2 x 10 -4 , as 
Mihos & Hernquist suggest that realistic values lie around 

10" 4 and less than 10" 3 . 



2. 1.2 Set-up for initial disc galaxies 

The baryonic mass of our model galaxies is 10 11 Mq. In dry 
merger simulations, this mass is purely stellar. In wet merger sim- 
ulations, 80% of this mass is stellar, and 20% is gaseous. The ini- 
tial gas and stellar discs are Toomre discs, with a scale-length of 
4 kpc and a truncation radius of 10 kpc for the stars, resp. 8 kpc and 
20 kpc for the gas. 20% of the stars are in a spherical bulge, mod- 
eled with a Hernquist ( 1990) profile with a 700 pc sc ale-length. Th e 
dark matter halo is modelled with a Burkert profile dBurkert|[l995l) . 
a 7 kpc scale-length and a truncation radius of 70 kpc, inside which 
the dark matter mass is 3 x 10 11 Mq. 

The two "progenitor" disc galaxies in each simulation are 
identical, the total mass of the remnant will be 2 x 10 11 Mq 
which is consistent with the slow-rotator mass range observed in 
the ATLAS 30 sample temsellem et al.ll2007l) . 

2.1.3 Orbits 

We have used two interacting orbits, for each kind of merger (dry 
and wet) and each resolution level. None corresponds to a very spe- 
cific and unlikely configuration like coplanar discs, or polar orbits. 

The first orbit is called "fast" because it forms fast- 
rotating early-type galaxies. The velocity at an infinite distance is 
170 km s _1 and the pericentre distance is 30 kpc. This orbit is pro- 
grade with respect to the first progenitor disc, with an inclination of 
the orbital plane w.r.t. the disc plane of 25°. The orbit is retrograde 
w.r.t. the other progenitor disc, with an inclination of 45°. 

The second orbit is called "slow" because it forms slow- 
rotating early-type galaxies (at least at high-enough resolution). 
The velocity at an infinite distance is 140 km s _1 and the pericen- 
tre distance is 25 kpc. This orbit is prograde with respect to the first 
progenitor disc, with an inclination of the orbital plane w.r.t. the 
disc plane of 45°. The orbit is retrograde w.r.t. the other progenitor 
disc, with an inclination of 25°. 

These orbits as well as those used in additional tests (Sec- 
tion 5) have a total energy E > or E ~ 0, correspond- 
ing to initially unbound galaxy pairs. Such orbits are repre- 
sentative of the most com mon mergers in A-CDM cosmology 
dKhochfar & Burkertll2006l) . 



Label 


Softening Length 


Particles / component 


Total particles 


standard 


180 pc 


10 b 


6 X 10 b 


high 


80 pc 


10 e 


6 X 10 e 


very high 


32 pc 


6 X 10 fa 


3.6 x 10 Y 



Table 1. Label for the resolution, softening length, number of particles per 
component (stars, gas and dark matter) and total number of particles in the 
simulation for the three resolutions. 
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Figure 1. Top panel \r profiles of all 200 projections for the Dry-Fast- 
very High simulation. The profiles are plotted as a function of R/R e (ra- 
dius normalised by the effective radius for each projection). Each line here 
corresponds to a given projection. Bottom panel Corresponding median 
(thick solid line), quartiles (thin solid line) and maximal and minimal values 
(dashed lines) at each radius. This representation of the results illustrates the 
fact that all values are between the dashed lines, and 50% of the projections 
are in the filled area. Each line plotted in this panel does not correspond 
to one specific projection: the median or quartiles are derived for different 
projections at each radius. 

2.1.4 Standard-, high- and very high-resolutions 

Wet and dry mergers have been simulated for each orbit at three 
resolution levels. The detail for these resolutions are indicated in 
TableQ] The very high resolution arguably corresponds to the high- 
est resolution simulation of a wet maj or merger performed so far 
(see lBournaud. Due & Emselleml2008h . 

We will label each simulation with the following nomencla- 
ture: 
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• the first item indicates a Wet or Dry merger, i.e. gas-rich or 
gas-free progenitors. 

• the second item specifies the chosen orbit: the one producing 
fast-rotators or slow-rotators (at least at high-enough resolution). 

• the last item indicates the resolution level: standard, high or 
very high. 

For instance, the wet-fast-high simulation refers to the high- 
resolution models of a wet merger on the orbit producing a fast- 
rotating early-type galaxy. 

2.2 Analysis of the relaxed merger remnants 

We analyse the merger remnants after 1.2 Gyr in the simulation, 
which is 800-900 Myr after the first pericentre passage, and 600- 
700 Myr after the central coalescence. The remnants are thus re- 
laxed when the analysis is performed. Tidal debris can still be or- 
biting around the merger remnant, but the bulk of the stellar mass 
in the central body does not show significant evolution. Analysis 
performed at earlier and later instants did not show significant vari- 
ations, so spurious effects related to time evolution should not affect 
the comparison of the three different resolution levels. 

2.2.1 Projected maps 

Intrinsic and apparent properties of the merger remnant (e.g the 
apparent ellipticity) are direct ly linked with its orbital structure 
( Jesseit, Naab & Burkert 2005). To probe the relaxed merger rem- 
nants, we have therefore built projected maps of the stellar mass 
density, line-of-sight velocity and velocity dispersion fields. Two- 
dimensional maps are useful to reveal the wealth of photometric or 
kinematic structures associated with a gala xy merger, e.g., globu- 
lar clusters or kine matic misalignments (see lBendo & Barne s2000; 
Ijesseit et al]|2007h . 

The projected maps cover a 16 x 16 kpc 2 field of view around 
the density peak of each system: our analysis is conducted up to a 
limit of three effective radii R e , which encloses most of the bary- 
onic mass of early-type galaxies, and the typical effective radius of 
our merger remnants is 2.5 kpc. Each projection was computed on 
a 100 x 100 pixel grid. The pixel size is 160x 160 pc 2 , which ap- 
proximatively corresponds to the size of the softening length of our 
standard-res simulations, and is kept fixed for all resolutions. 

To obtain statistically significant results, we have built such 
maps and performed the subsequent analysis with 200 isotropically 
distributed viewing angles (i.e., 200 different line-of-sights). In this 
way, we do not characterize and compare the merger remnant under 
a particular projection, but their global, statistical properties. As an 
example, Figure[T]shows the effect of the projections on the radial 
\r profiles for one simulation. Among these 200 profiles, the low- 
est (near zero) and highest values correspond respectively to the 
merger remnant seen nearly face-on (i.e., the lowest apparent ellip- 
ticity) or nearly edge-on (i.e., the highest apparent ellipticity). Our 
choice of 200 projections ensures that neighbouring projections are 
separated only by 10 degrees in any direction, so that intermediate 
viewing angles would not show significant differences. 

2.2.2 Physical parameters 

Our analysis is based on a few simple morphological and kinematic 
parameters - a choice mainly motivated by the fact that these pa- 
rameters are often being used as standards in studies of nearby el- 
liptical galaxies. 



The morphological parameters pertains to the photometry : we 
measure the ellipticity e (defined as 1 — b/a, where a and b are 
the semi major- and minor-axes, respectively) and a^/a which is 
the fourth (cosine) Fourier coefficient of the deviation of isophotes 
from a perfect ellipse (a^/a > and a^/a < correspond to discy, 
and boxy isophotes, respectively). These two parameters are com- 
puted using the Kinemetry software tooQ which can be used to per- 
form standard elli pse-fitting of galaxy im ages, as well as to study 
galaxy kinematics IKrainov ic et alj2006l) . For the kinematic analy- 
sis, apart from the first two velocity moments (velocity and velocity 
dispersion), we use the Ah parameter, a robust pro xy for the bary- 
onic p rojected angular momentum, as defined in Emsellem et al. 
j200l hereafter E07) : 

A _ (R\v\) 

In E07, Xr was used to reveal two families of early-type 
galaxies, the slow-rotators with Ah ^ 0.1 and the fast-rotators 
with Xr > 0.1 a t one effective radius R e . In a recent study, 
Ijesseit et al.l J2008h have simulated binary disc mergers to inves- 
tigate the Ah parameter: tests on their merger remnants reveal that 
An is a good indicator of the true angular momentum content in 
earl y-type galaxies. As em phasised in E07. ICappellari et all J2007h 
and IKrainovic et al ] d2008l) . fast and slow rotators exhibit qualita- 
tively and quantitatively different stellar kinematics. Ah is thus an 
interesting parameter to probe, and should indicate whether or not 
the kinematics of the merger remnants are equally resolved at dif- 
ferent resolutions. 

For each above-mentioned parameter, we have computed the 
minimum, maximum, mean values, as well as the I st and 3 rd quar- 
tiles over all the projections at individual radii, to quantify the sta- 
tistical distribution of these parameters in a simple way. An exam- 
ple is shown in FigureQ] Note that with this choice, the projection 
which minimises or maximises a parameter varies with radius. 



3 EFFECT OF RESOLUTION ON THE FORMATION OF 
SLOW ROTATORS 

In this Section, we briefly describe the properties of the simulated 
mergers with the three different resolutions. The complete set of 
analysis results can be found in Appendix A. We then focus the 
analysis on the simulations that show important differences, namely 
the cases producing slow rotators at high resolution. 

3.1 Morphology and kinematics 

Radial stellar density profiles are shown in Figure|2] We then show 
in Figure[3]the projected stellar density maps, of the relaxed merger 
remnants in all simulations, choosing the flattest and roundest pro- 
jections as well as a projection representative of the mean ellipticity 
in each case. The corresponding line-of-sight stellar velocity fields 
are presented in Figure [4] f or the same projectio ns, the maps have 
been Voronoi binned (Cappellari & Copin 2003) to the same level 
of 15 particles minimum per bin. Further morphological or kine- 
matics parameters are presented in Appendix A. 

This analysis reveals various similarities or differences, de- 
pending on which merger is considered. The most noticeable re- 
sults are as follows: 

1 http://www-astro.physics.ox.ac.uk/ dxk/idl/ 
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Dry-Slow 

Emax tmeart 



standard-res 



high-res 



very high- res 




Wet-Slow 




standard-res 



high-res 



very high-res 



Figure 3. The twelve normalized projected surface density maps (in log), for the four sets of simulations at three different resolutions (labelled accordingly). 
The field of view is 16 X 16 kpc . For each simulation, the projections corresponding to the minimum, maximum and mean ellipticities are shown. The 
viewing angle of these projections are defined at very high-res and re-applied for the standard-res and high-res simulations: projections are thus established 
along the same line-of-sights for all resolutions. Luminosity contours are the same for all simulations and drawn with a spacing of 0.5 magnitude (except for 
the 2 inner contours with a step of 0.3). The effective radius is about at the edge of the 4 th isophote for all simulations. 



• Mergers that produce fast-rotators at the highest resolution 
also result in fast rotating systems at the lower, standard resolution. 
Overall, the apparent morphology for any projection of the Dry- 
Fast and Wet-Fast models is unaffected by the resolution (Figure[3]l. 
The velocity fields are also quite similar (Figure[4]l, with only mi- 
nor misalignments between the apparent kinematic and photomet- 
ric axes. Ellipticity and \r profiles, provided in Appendix A (see 



Figure IA1I and IA2t . confirm these similarities and that all these 
mergers remnants are fast rotators, with a rotational support that is 
largely independent from the numerical resolution. 

• Strong kinematic misalignments and kinematically decoupled 
cores (KDCs) are found only in slow-rotators, but really appear 
only at high resolution. The Dry-Slow model has a KDC at stan- 
dard resolution, but its amplitude is significantly lower than the one 
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observed in the high and highest resolution models. The Wet-Slow 
model has a KDC only at high/very high resolution. Overall, kine- 
matic misalignments increase at high resolution, as illustrated for 
instance by the flattest projections of the Wet-Slow case. 

• Morphological and kinematic differences are most important 
for mergers that produce slow -rotators at high resolution. Striking 
morphological differences are seen in particular for the Wet-Slow 
case (Figures|2]and[3j and both the amplitude and the shape of the 
velocity field change with resolution for the Wet-Slow and Dry- 
Slow cases (Figure |4j- For instance a rapidly rotating core is seen 
in the Wet-Slow merger remnant at standard-resolution, instead of 
a slow-rotating KDC at high- and very high-resolutions. The Dry- 
Slow remnant also shows up as a discy rotating system at standard 
resolution, in contrast with the observed remnant at higher resolu- 



tions. We also note on Figure [2] that the stellar density profile is 
resolution-dependent in particular for the Wet-Slow case, with a 
much less concentrated merger remnant in the standard-resolution 
case (the mass within 5 kpc is about 25% lower than at high or very 
high resolutions). 



3.2 Formation of slow-rotators at high resolution 

We now focus on the detailed properties of the mergers for which 
the most important differences have been noticed, namely those 
producing slow rotators at the highest resolutions. 
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R (kpc) 




Figure 2. Stellar density profiles for the Wet-Slow, Wet-Fast, Dry-Slow and 
Dry-Fast simulations (from top to bottom, respectively). The density of the 
Wet-Fast, Dry-Slow and Dry-Fast cases has been divided by a factor of five, 
ten and twenty respectively to improve the readability of the plot. Red lines 
correspond to very high-resolution models, green lines to high-resolution 
cases, and blue lines to standard-resolution models. The two alternative re- 
alisations of the Wet-Slow-Standard simulation (see Sect. 14.31 are shown in 
dashed and dotted lines. 



Old stars 



Young stars 



Gas 




Figure 5. Projected density maps of the old stars, young stars, and gas in 
the wet-slow remnants; from top to bottom : standard-, high- and very high- 
res models, for the projection which minimises the ellipticity (e m i n ) as in 
Figure [3] Old stars are those formed before the merger, young stars are 
formed during/after the merger. The field of view is 16 X 16 kpc 2 , and the 
isocontours correspond to the projected old stellar component. 




3.0 



Figure 6. Top panel Radial line-of-sight velocity and velocity dispersion 
profiles (respectively bottom and top lines of the plot) for the mean ellip- 
ticity projection along the global kinematic position angle of the Wet-Slow 
simulation (right panels of Figur^4). Bottom panel A^ profiles as a func- 
tion of R/ R e , the minimum, median, maximum and quartiles values are 
presented as in Figure [T] In both panels, the standard-res is represented in 
blue, the high-res in green, the very high-res in red. 




Figure 7. A^ profiles as a function of R/R, e for the Dry-Slow simulation. 
The standard-res is represented in blue, the high-res in green, the very high- 
res in red. 
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3.2.1 Morphology and Kinematics 

To better understand the differences seen in the morphology of the 
Wet-Slow simulations, we have examined the three included bary- 
onic components of the merger remnants separately, namely the 
"old" stars formed before the beginning of the merger event, the 
"young" stars formed during/after the merger event, and the gas 
left over after the merger (see Figure[5j. Within the central 10 kpc, 
the standard-res remnant exhibits a prominent bar, the inner distri- 
bution of the gas and young stars being driven by this tumbling 
structure with e.g., a ring-like structure at a radius of ~ 6 kpc. 
In the high-res and very high-res, the gaseous component and the 
young stars have a smoother distribution more closely following the 
overall old stellar distribution. In addition, many young star clus- 
ters are visible in the maps from the very high-res, a few in the 
high-res, and none in the standard-res. High spatial resolution of 
course allows to resolve the forma tion of stellar clusters (see also 
iBournaud. Due & Emsellemll2008l) but there is also a larger num- 
ber of other young stellar substructures at increasing resolution, 
like filaments, tidal streams, and a compact nucleus (Figure [5j. 

The kinematic discrepancies discussed above in the velocity 
fields are quantified globally in the radial velocity and Xr profiles 
(Figure [6). The standard-res displays significant rotation inside 3- 
4 kpc (up to ~ 85 km s _1 ) and a decreasing rotation velocity at 
larger radii. There is a drop in the velocity dispersion in the central 
2 kpc, and no sign of a KDC. This is in stark contrast with both 
the high-res and very high-res which overall show much lower ro- 
tational velocity support (below ~ 50 km s _1 and particularly low 
in the central 2 kpc), and a KDC in the central 1 kpc. Overall the 
high-res and very high-res have similar velocity rotation curves, 
apart from a more pronounced KDC signature in the very high-res 
(partly due to the KDC having a slightly different position angle in 
these two remnants). 

The general discrepancies of the standard-res versus high-res 
and very high-res realisations are confirmed by the Ah profiles 
(Figure [6](. The merger remnant made at standard-res is clearly a 
fast-rotator. The high-res and very high-res are both classified as 
slow-rotators with respectively a maximum value of A_r of 0. 1 and 
0.06 at one effective radius. The Xr profile goes up somewhat more 
rapidly with radius in the high-res case than in the very high-res, but 
the difference remains of the order of the scatter between different 
projections of each case. The presence of a bar in the stellar compo- 
nent of the standard-res is likely a result of the significantly higher 
rotational support (see also SectO. 

Beyond one R e , the Xr profiles of the high-res and very high- 
res are rising: there is less angular momentum in the center, which 
has been expelled outwards (see also E07). However, even at these 
large radii, the slow-rotators have less angular momentum than 
fast-rotat ors (see FigureRp. Observations conducted up to two or 
three R e JWeiimans et alT2009l ; ICoccato et ai]|2009t) would bring 
additional constraints on the formation scenario of slow-rotating 
early-type galaxies. 



3.2.2 Role of gas on the properties of merger remnants 

The Dry-Slow simulations show smaller differences in the stellar 
density maps and velocity fields. They also exhibit smaller differ- 
ences in their Xr profiles (Figure[7}. Nevertheless, the standard-res 
simulation is again a faster-rotator than the high-res and very high- 
res cases at 1, 2 or 3 effective radii. A KDC is also found only in the 
high-res and very high-res cases, associated to a peak of Xr inside 
one effective radius. 



A lower specific angular momentum in the main stellar body 
at higher resolution is not only found in Wet-Slow mergers, but 
also in Dry-slow mergers, the differences being still much more 
pronounced in the Wet case. 

Ga s plays an important role i n shaping merge r rem- 
nants jNaab Jesseit & Burkertl 120061 : iRobertson et al.l l200d 
iHopkins et aU l2009b ) and it is interesting to compare the Wet- 
Slow and Dry-Slow merger remnants at fixed resolution, to better 
understand its specific impact (Figures[6]and[7J: 

• at standard-res, the wet merger remnant has a much higher 
rotational support than the dry case. This is consistent with the 
usually known effect of gas helping the survival of rotating stel- 
lar discs during major mergers, and/or re-building of discs after 
mergers (Robertson et al. 2006, Hopkins 2009b). 

• at high-res and very high-res, the rotational support of the 
merger remnant is not increased when gas is present. The angular 
momentum, traced by Xr, is actually lower by about 20% inside 
one effective radius in the very high-res wet case, compared to the 
corresponding dry merger. 

It thus seems that the impact of gas on the global properties of 
major merger remnants is more complex than originally thought, 
and can even be weakened at high resolution. This suggests that 
the global dynamics of gas during the major merger or in a young 
merger remnant can be significantly affected by resolution. As seen 
in Figure|5] gas at standard-res largely lies in smooth structures and 
the formation of new stars during the merger proceeds in a rela- 
tively smooth way. At increased resolutions, thinner gas structures 
are resolved during the merger, which can result in clustered star 
formation and the formation of numerous young stellar structures, 
as observed in the final merger remnant in Figure|5] 



3.3 Summary of the resolution tests 

The resolution does not seem to significantly affect the morphol- 
ogy and kinematics of the mergers remnants that are fast rotators at 
high resolution: they are still fast rotators at lower resolution, with 
very similar morphological and kinematic properties. This con- 
trasts with the fact that resolution has a major effect on the forma- 
tion of slow-rotating systems. The systems that are slow rotators at 
high resolution rotate more rapidly when the resolution decreases, 
and can be observed as true fast rotators at standard-res. The effect 
is small in dry mergers, but is dramatic in our wet merger model. 
KDCs in these slowly-rotating systems are also significantly bet- 
ter resolved at high resolution. The role of gas in shaping merger 
remnants is found to vary with resolution: at low resolution, gas 
rebuilds rotating disc components, increasing the overall disciness 
and rotational support. At higher resolution, the effect cancels out: a 
merger that forms a slowly-rotating system in a dry case still forms 
an equally slow-, or even a bit slower-rotator in the corresponding 
wet case. 

The next section focusses on interpreting the origin of the res- 
olution effect in the formation of slowly-rotating ellipticals. We in 
particular show that it is not an artifact caused by different initial 
conditions or a bias in the simulated orbits, but a real effect related 
to the way the violent relaxation during the merger itself is treated. 



4 ORIGIN OF THE RESOLUTION EFFECT 

We here show that the above-mentioned discrepancies observed in 
the simulations that produce slow-rotating ellipticals at high resolu- 
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Figure 8. Xr profiles of the two Wet progenitors as a function of R/R e . 
Colours, as in previous Figures, with the standard-res in blue, the high-res 
in green, the very high-res in red. 

tion, are really attributable to the physical modeling of the merging 
process. They are not artefacts related to initial conditions of the 
progenitor galaxies and/or interaction orbits that would vary with 
the resolution. 



4.1 The progenitor galaxies 

We first check that the progenitors galaxies are similar at any res- 
olution. To this aim, we analyze their kinematic properties, in par- 
ticular the Xr profiles - e and a^/a parameters are less relevant for 
disc-dominated galaxies. Since the merger simulations were per- 
formed after an isolated relaxation of each progenitor galaxy (see 
Sect. |2j, we analyzed the progenitors from a snapshot right after 
this relaxation period, so that the results (Figure [8} are representa- 
tive of the conditions under which the mergers occur. 

The two progenitor galaxies have quite similar angular mo- 
mentum profiles (Figure[8}. There are some fluctuations, but they 
are not systematically corresponding to an increase or decrease of 
Ah with resolution. They are also weaker than the discrepancies 
found in the final merger remnants. Actually, they result for a large 
part from the effective radius changing slightly with the resolution, 
and profiles of Xr as a function of the absolute radius (in kpc) show 
smaller differences than the profiles in units of the effective radius. 
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Figure 9. e (top) and Xr (bottom) profiles in function of R,/R e for 5 
slightly different orbits for the Wet-Slow-High simulation. 

These fluctuations cannot therefore be the main cause for the ob- 
served resolution effects in the merger remnants. 

4.2 Interaction orbits 

Simulations at the three resolutions are started with the same rel- 
ative position, velocity and inclination for the two interacting pro- 
genitors. However, varying the resolution may result in slight dif- 
ferences in dynamical friction and angular momentum exchanges, 
if these processes are resolved differently, and the interaction orbits 
might diverge before the merger actually takes place. If this were 
the case, our results would be attributable to different orbits rather 
than different treatments of the merging process itself. 

We found that the positions at the first pericentric passage 
vary by 2.1 kpc on average, and the velocities by 9 km s _1 . Al- 
though these differences seem small and no systematic variation 
with resolution appeared, we further investigated their potential ef- 
fect. To this aim we performed four new realizations for the Wet- 
Slow model at high-res, with variations of the position or the ve- 
locity twice larger than the average values above (i.e., ±3.6 kpc 
and ±18 km s _1 , respectively). The results are shown in Figure [9] 
for the morphological and kinematic profiles of e and Xr. Changes 
are minor and differences arising in the interaction orbits cannot 
explain the variation of the results with resolution. 
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Figure 10. Top panel Stellar (old plus young) intensity maps for three dif- 
ferent realisations of the Wet-Slow-Standard simulation. The field of view is 
16 X 16 kpc 2 . The simulation used in the study is in the left panel. Bottom 
panel A^j profiles as a function of R/R e of the three above simulations. 



4.3 Robustness of the Wet-Slow-Standard simulation 

As the resolution effect found in the formation of slow rotators, in 
particular in the Wet-Slow model, cannot be attributed to a change 
in the initial conditions and interaction orbit, it likely relates to the 
physical treatment of the merging process itself. Nevertheless we 
wanted to check whether or not this could still be attributed to the 
presence of particle noise, which is higher in the standard-res cases. 

The Wet-Slow model at the standard-res shows a strong stel- 
lar bar, contrary to the high-res and very high-res cases. We wanted 
to make sure that this bar is a robust consequence of the high rota- 
tional support of the standard-resolution case, and is not a misinter- 
preted effect that arised from a particular realization of the particle 
noise. 

To this aim, we performed two other Wet-Slow-Standard sim- 
ulations with the same initial conditions but different, random re- 
alizations of the particle noise. The final stellar distribution, shown 
in Figure [10] all show a similar bar, and the profiles are also 
relatively similar to the original Wet-Slow-Standard model - there 
are some variations, but the A_r distributions of the three realiza- 
tions overlap with each other, and the three systems are equally fast 
rotators. These two new realisations are also shown in dashed and 
dotted lines on Figure l2l and 1111 and again share common proper- 
ties with the initial Wet-Slow-Standard model, and hence the same 
differences compared to the higher-resolution cases. 

Thus, the role of bars and spiral patterns in redistributing the 
mass and angular momentum in the standard-resolution Wet-Slow 
model is robust, independent of a particular realization of the parti- 
cle noise. We also find (see next subsection) that the time variations 
of the gravitational potential during the interaction and merger are 
similar for the three standard-res realizations. 
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Figure 11. Top panel Median variations of the gravitational potential (in 
arbitrary units) over 5000 test particles as a function of time. The two other 
realisations of the Wet-Slow-Standard simulation are shown in dashed line. 
Top panel Histogram of the maximum variation of the gravitational poten- 
tial (in arbitrary units) of each particle over the simulation. 



4.4 The on-going merger 

At this point, we have established that the differences observed 
in the merger remnants do not result from variations in the ini- 
tial conditions, interaction orbits, or particle noise. The differences 
should then arise in the physical treatment of the merging process, 
which would mean that they are "robust" effects, potentially af- 
fecting any simulation with any numerical code. Varying the spa- 
tial and mass resolution could affect the detailed evolution of the 
dissipative component (including star forming structures), and this 
could in turn modify the overall orbital structure of the merger rem- 
nant (seelBarnes & Hernguistlll996|; iNaab, Jesseit & BurkerfeoOfj; 
ICox et al.ll200q) . However we have seen that the resolution effect 
does not completely disappear in dry mergers. A more general ef- 
fect can be the treatment of the violent relaxation, i.e. the rapid 
changes of gravitational potential that are responsible for the evac- 
uation of energy and angular momentum from the main body of 
the merger remnant - these quantities being carried away by a low 
fraction of the mass expulsed at large radii. This process of course 
plays a more important role in the formation of slow rotators than in 
the formation of fast rotators. The resolution effects are much more 



© ?? RAS, MNRAS Q00.nHT4l 



Formation of slowly rotating early-type galaxies 1 1 



important for slow rotators than fast ones (Sect. [3}, which suggests 
that they do actually relate to the violent relaxation process. 

To quantify the importance of violent relaxation in our merger 
simulations, we followed, in the Wet-Slow models, the variations of 
the gravitational potential D<j>/ Dt of 5000 randomly chosen "test" 
particles, all of which are stellar particles existing at t = 0, all 
along the simulations. The derivate is Lagrangian, since it follows 
the motion of each particle. In an isolated galaxy, D<j>/ Dt relates to 
the variation of potential along the orbit of each particle, in partic- 
ular their radial excursion in the potential well of the galaxy. Dur- 
ing the interaction and mergers, peaks of D<j>/ Dt should trace the 
importance of scattering by local density fluctuations through the 
violent relaxation process. 

The top panel of Figure [TT] shows the median value of 
\D<f>/Dt\ as a function of time - we take the absolute value for 
each particle, as a particle moving inward or outward can be con- 
sidered with the same behaviour. Before the merging (i.e before 
~150 Myr) the three resolutions are identical, meaning that there is 
no difference in the progenitors during the approach phase, modest 
values of \D(f>/Dt\ simply correspond to modest radial excursions 
of particle in the progenitor disc galaxies. 

After the merger, each simulation shows a relatively constant 
| D<f> I ' Dt\ in a relaxed system, but the value is higher at high-res and 
very high-res, indicating larger radial excursions of stellar particles 
compared to the standard-res case. More radial orbits are indeed 
expected for slow rotators compared to the standard-res fast rota- 
tor. We note again that the different orbital structure does not only 
affect the gas and the young stars formed during the mergers, but 
also the old stars present before the merger itself (see also Figure[5] 
and Sect.[3j. 

During the merging process, a first peak in the median 
\D(f>/Dt\ occurs at the first pericenter passage, after about 
150 Myr, but is more pronounced at high(est)-resolution. Another 
peak is found at the high-res and very high-res during the final co- 
alescence at t ~ 280 Myr , but is much weaker in the standard-res 
case. The final coalescence does take place at the same moment for 
the three resolutions, but is a smooth process in the standard-res 
case, while it is accompanied by rapid variations of the potential 
undergone by stellar particles at high resolution. The bottom panel 
of Figure [TT1 shows the maximum variation of \D(f>/Dt\ for each 
particle over all the simulation. The distribution at standard-res is 
very different from the distribution at high-res and very high-res. 
This confirms that the particles at standard-res undergo less rapid 
variation of the potential, i.e. lower peaks of gravitational forces. 

This overall demonstrates that the relaxation process, during 
the merging of galaxies, is smoother at low resolution than at high 
resolution. We have shown previously that the high-res and very 
high-res simulations resolve much more dense substructures, like 
gas filaments, stellar clusters, compact cores, etc. Our interpreta- 
tion is then that these local density peaks are accompanied by rapid 
variations of the gravitational potential, which scatter the stellar 
orbits, evacuate the angular momentum, and form, for favorable 
orbits, slowly-rotating elliptical galaxies. At low resolution, these 
rapid and local fluctuations of the density and potential are largely 
missed, hence the merging process is smoother, and more angular 
momentum remains in the main stellar body of the merger remnant. 

Density fluctuations are of course stronger in the dissipative 
component (gas) and the young stars formed therein, which likely 
explains why the resolution effect is stronger in wet mergers. Nev- 
ertheless, old stars are clearly affected as well, as was shown above. 

This also explains why the effect of gas in a wet merger, 
compared to a dry merger at fixed resolution, is different for 



standard-resolution models and high-resolution ones (Sect. PO^I l. 
At standard-res, the gas remains relatively smooth, promotes the 
survival/rebuilding of a stellar disc component, thus increasing the 
rotational support in the final merger remnants. At higher reso- 
lution, the presence of gas forms many dense small-scale sub- 
structur es of gas and young stars (consi stent with observations, 
see e.g. lBournaud. Due & E msellem 2008), these substructures in- 
crease the degree of relaxation during the merging process, not just 
for the gas and young stars but also for the old stars. Thus, while 
the presence of gas should still promote the survival/rebuilding of 
a disc component in the merger remnant (our high-resolution wet- 
slow remnant does have a low-mass disc componnent of gas and 
young star), it also promotes orbital scattering and evacuation of 
the angular momentum for the whole baryonic mass, but the latter 
effect is missed if the resolution is too low. This explains why, at 
high resolution, the merger remnant (in the Wet-Slow case) does 
not have a higher rotation support or a more prominent disc com- 
ponent than the corresponding Dry-Slow case, an in fact even has a 
somewhat lower \r at one effective radius. 

The high-resolution simulations, compared to the standard 
cases, resolve the formation of dense and relatively massive sub- 
structures (clusters, cores, filaments of 10 5 ~ 7 solar masses) that 
scatter the stellar orbits and evacuate the angular momentum from 
the main body of the merger remnant. Very high-res simulations 
show a relatively reasonable convergence compared to the high-res 
ones: they resolve the same massive substructures, plus lower mass 
ones (~ 10 4-5 solar masses) that are more numerous but are much 
less efficient to scatter the orbits and affect the relaxation of the 
merger remnant, as the corresponding relaxation timescale is much 
longer. It is thus expected that results converge at a high-enough 
resolution. 

4.5 Timestepping and code specificities 

Our results have been obtained with a given code and one can nat- 
urally wonder whether or not other codes would show the same 
resolution effect. We in fact expect no fundamental differences in 
the output from different codes, given that similar substructures are 
formed initially: this relaxation effect is mostly gravitational, and 
this should be treated rather similarly in grid-based and tree-codes. 
The main question remains then whether or not other codes would 
form substructures similar to those found in our simulations (e.g ., 
with a similar mass spectrum. iBournaud, Due & Emse llem 2008): 
this a priori depends on the modeling of gas cooling and turbulence 
dissipation processes. 

Another specificity of the code employed is its fixed timestep. 
A small timestep may better resolve the scattering of stellar or- 
bits by dense substructures, in particular at high resolution. This 
would actually tend to increase the effect of resolution that we 
have found, which justifies studying the resolution effect at fixed 
time-step rather than decreasing the timestep at increasing spatial 
resolution. This way, the effects found can be robustly attributed 
to the spatial resolution. The time-step itself may have additional, 
separate effects, in our code or any other, that should be studied 
separately at fixed (high) resolution. 

5 A SYSTEMATIC EFFECT IN THE FORMATION OF 
SLOW ROTATORS 

To ensure that the resolution effect in the formation of slow rotators 
is a systematic one, and not specific to one simulated orbit, we have 
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selected other major mergers that form slow rotators at high-enough 
resolution in a larger simulation sample (Bois et al. in preparation), 
and re-simulated them at lower resolution. These three additional 
mergers were not simulated at the very high-res but at a resolution 
which is actually a bit higher than the high-res, with a spatial res- 
olution of 58 pc and 2 x 10 6 particles per component par galaxy 
(i.e. a total of 1.2 x 10 7 particles). They were also resimulated at 
the same standard-res level as the previous models. The gas frac- 
tion in these models is 10%, so as to check that a high gas fraction 
is not required for the resolution effects to arise. Compared to the 
original slow-rotator simulation, one orbital parameter is changed 
in each case. Model slow-1 has a retrograde orbit for both galaxies, 
model slow-2 has a pericenter distance of 25 kpc, model slow-3 has 
an orbit inclination of 25 degrees for both discs. 

The projected stellar density maps and line-of-sight veloc- 
ity fields are shown for these models, at standard- and high- 
resolutions, under the flattest projections in each case, on Fig- 
ure [T2] (the Xr profiles being also shown in this figure). These 
three merger remnants are slow rotators (at least at one effective 
radius) at high-res, but at standard-res they all have a much higher 
angular momentum \r, and a velocity field aligned with their mor- 
phological axis - the high-resolution slow rotators have important 
kinematic misalignments and KDCs. 

These additional cases confirm the effect found and analyzed 
in detail in the original slow-rotator model, with about the same 
degree of discrepancy between the standard-resolution and higher- 
resolution models. Because we have spanned the four slowest ro- 
tators among the 1:1 mergers from the Bois et al. (in preparation) 
sample, the resolution effect seems to strongly affect the modeling 
of a significant number of slow-rotating ellipticals, if not all, even 
with modest gas fractions (here 10% of the baryonic mass). 



6 DISCUSSION AND CONCLUSION 

In this paper, we have studied the effect of numerical resolutions 
(spatial and mass resolution) on the global properties of merger 
remnants. Our simulations at "standard" resolution, are comparable 
to the majority of merger simulation samples published in the last 
few years: the spatial resolution (gravitational softening and typi- 
cal hydrodynamical smoothing lengths) is 180 pc, and the number 
of particles ~ 10 5 per galaxy and per component (gas, stars, and 
dark halo). These simulations have been compared to models of 
the same mergers with increased resolution, up to 32 pc and almost 
10 7 particle per galaxy and per component. 

We have analyzed the morphology and kinematics of the re- 
laxed merger remnants. In particular, we have studied whether they 
are "fast rotators", with an apparent spin parameter Xr > 0.1 and 
have small misalignments between the morphological and kine- 
matic axes, i.e. in broad terms early-type galaxies with significant 
flattening and rotational support. At the opposite end, "slow rota- 
tors" are systems with a low Xr ^0.1 (at one effective radius), 
large kinematic misalignments, i.e. early-type galaxies dominated 
by (anisotropic) pressure support and low residual rotation. Such 
slow rotators usually have central KDCs in our high-enough reso- 
lution simulations. 

Our mains findings can be summarized as follows: 

• The formation of fast-rotators is not significantly affected by 
numerical resolution. Models that produce fast rotators at the high- 
est resolution also result in fast rotators at lower resolution, with 
some random fluctuation of their properties but no sign of sys- 
tematic variation in the morphology or angular momentum profile 
against resolution. 

• The formation of slow-rotators is greatly affected by numeri- 
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cal resolution. Models that produce slow rotators at the highest res- 
olution result in much faster rotators at lower, standard resolution. 
The effect is present, but relatively minor, in purely collisionless 
dry mergers. Discrepancies become major in wet mergers, even in 
cases with modest gas fractions like 10% of the baryonic mass. 

• These effects cannot be attributed to our choice of initial con- 
ditions or interaction orbits, but actually relate to the physical treat- 
ment of the merging process itself. In particular, small-scale den- 
sity fluctuations increase at high resolution, and they participate to 
scattering stellar orbits and largely influence the final degree of re- 
laxation and orbital structure in the merger remnants. 

• The effect of gas on the properties of merger remnants is gen- 
erally considered to consist in preserving a higher angular momen- 
tum, in particular through enhancing the survival/rebuilding of disc 
components in merger remnants. We find that this picture is incom- 
plete: at high resolution, gas still reforms discy components, but 
also forms a large number of dense substructures (massive star clus- 
ters, dense nuclei, tails and filaments, etc) that trigger rapid varia- 
tions of the gravitational potential and the degree of relaxation of 
the final system. This effect is missed with a too low resolution. 
At high-enough resolution, adding gas to a given merger does not 
necessarily increase the rotational support of the final merger rem- 
nant; we even find a case of a wet merger with 20% of gas that has 
a final angular momentum parameter Ar slightly lower than the 
corresponding dry merger. 

At the present stage, our results do not indicate how frequently 
real slowly-rotating ellipticals were formed by binary (wet) merg- 
ers of disc galaxies, but they show that this can be a robust path- 
way for their formation. In the course of the ATLAS^ D project, we 
are conducting, analyzing, and comparing a large set of numerical 
simulations for various formation mechanisms, in order to derive 
which is (are) the main formation mode(s) of real slow rotators in 
the nearby Universe. Our present results already indicate the limi- 
tations of existing samples of galaxy merger simulations, and will 
then serve to estimate the required resolution, the limitations of nu- 
merical models and their possible biases. 

More generally, the immediate implications of these findings 
on our understanding of early-type galaxy formation are: 

• High resolution in simulations of major mergers does not just 
allow to resolve small-scale structures like nuclear systems and star 
clusters, but impacts the whole global properties of the elliptical- 
like merger remnants, at least for the slow-rotating ones. 

• The formation of slow-rotating elliptical galaxies can be 
achieved through a major merger relatively more easily than pre- 
viously believed. It can be frequent even in wet mergers with rela- 
tively high gas fraction, and with late-type, disc-dominated progen- 
itor galaxies. 

• Repeated mergers and/or dry mergers of galaxies that are al- 
ready early-type systems are thus not the only theoretical path to 
produce slow-rotating galaxies. Major mergers of two die galaxies, 
including wet mergers, can produce slow-rotating early-type galax- 
ies. Further studies are needed to determine how common this for- 
mation mechanism is for slow-rotators. 

• Quantitative comparisons of major merger simulation results 
with the observed properties of real early-type galaxies require high 
resolution models. A typical requirement, according to our study, 
would be a spatial resolution better than 100 pc for both the grav- 
itational N-body aspects (i.e., softening length) and the hydrody- 
namical ones (for instance, the size of groups of particles other 
with quantities are smoothed in SPH models). The mass resolution 
should correspond to at least ~ 10 6 particle per galaxy per compo- 



nent, which typically corresponds to a mass resolution ~ 10 4 Mq 
for the gas discs of bright spiral galaxies. We find reasonable con- 
vergence above this resolution, but cannot rule out that some sys- 
tematic effects still exist; in any case simulations below this reso- 
lution level show clear and strong resolution effects. Unfortunately, 
most published samples of major merger simulations are below this 
typical resolution limit. 

• The small-scale properties of interstellar gas and clustered 
star formation are important for the global, large-scale proper- 
ties of merger remnants. Simulations directly resolving gas cool- 
ing down to low temperatures, the formation of cold (molecu- 
lar) gas clouds and star formation therein, are highly desirable 
to understand the whole process of early-type galaxy formation. 
Modern hydrodynamic codes are promising in this respect (e . g . 
iBournaud. Elmegreen & Martidfe oog; Kim. Wise & Abelll2009h . 
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APPENDIX A: ADDITIONAL MORPHOLOGICAL AND 
KINEMATIC PARAMETERS FOR THE MERGER 
REMNANTS 

In this appendix, we provide further details on the morphology and 
kinematics of the four merger remnants analysed in the present pa- 
per. 

Al Morphology 

The ellipticity e and 04 /a profiles are shown in Figure IaTI as a 
function of R/R e . The apparent differences sketched in Sect. 13.1 1 
are confirmed quantitatively in the radial ellipticity profiles. Within 
1.5 R e , there are small differences in the Dry/Wet-Fast and Dry- 
Slow simulations. The ellipticity outside 1.5 R e is however quite 
similar at all resolutions for these three simulations: the minimum 
ellipticity is basically 0, the mean is 0.33 ± 0.03 and the maximum 
is 0.55 ± 0.05 at 2 R e for all three standard-res, high-res, very 
high-res. 

The Wet-Slow simulation shows much larger differences. Out- 
side 0.5 R e , the high-res and very high-res are similar. The ellip- 
ticity profile of the standard-res has then a completely different 
appearance: between 0.6 and 2 R e , 75% of the projections have 
an ellipticity higher than 0.4, and the reached maximum in e is 
0.75 (versus ~ 0.6 for the other two resolutions). In the outer part 
(R > 2R e ) the ellipticity of most of the projections is decreasing 
but its maximum is still larger than 0.7. 

The same trends are observed in the 04/ a profiles. In the 
Dry/Wet-Fast and Dry-Slow simulations at all three resolutions, the 
mean 0,4 /a is around 0. Then 50% of the projections are between 
-2 and 2%. The high-res Wet-Fast and the very high-res Dry-Fast 
simulations are only slightly more boxy. The 0,4/0, profile of the 
Wet-Slow simulation dramatically confirm what we observe for the 
ellipticity. The standard-res clearly departs from the high-res and 
very high-res, which are quite similar. Between 0.5 and 1 Re, the 



projections of the standard-res span a very large range of 04/0. 
Between 1 and 2 R e , 75% of the projections have a discy shape, 
and the isophotes of the merger remnant become then increasingly 
boxy going outwards. 

A2 Kinematics 

The Wet-Slow simulation has been treated in the paper, we will thus 
focus on the three other simulations. 

Left panels of Figure lA2l show the velocity and velocity dis- 
persion curves for the mean-ellipticity projection along the global 
kinematic position angle. In the Dry-Fast simulation, the central 
slope of the rotation curve at very high-res is slightly shallower, and 
the dispersion about 15% smaller, but in the outer part the velocity 
amplitude is similar at all three resolutions, with a velocity of about 
60 km s _1 at 6 kpc, and dispersions values going to about 150 km 
s - . The Wet-Fast simulations show consistent velocity profiles at 
all three resolutions, with a velocity amplitude of 60 km s _1 at 
6 kpc, and dispersion decreasing outwards down to ~ 125 km s - . 
Again, the velocity curves for the Dry-Slow simulation are all very 
similar, but these profiles clearly reveal the previously observed 
KDC which appears here as a kpc-size core counter-rotating with 
respect to the outer part. Note the standard-res dispersion profile 
which is about 10% smaller than for the other two higher resolu- 
tions. 

In right panels of Figure lA2l we now compare the simulations 
using the apparent angular momentum Xr. These figures clearly 
show that the Dry/Wet-Fast simulations (top and second from top) 
both result in fast-rotators, the mean values of Ah is 0.2 and the 
maximum about 0.25 at 1 R e for the three resolutions. This con- 
firms our previously mentioned results that the spatial and mass 
resolutions do not seem to have a significant effect on these merger 
remnants. 

The analysis of the morphology and kinematics of the Dry- 
Slow simulation did show mild differences in the remnants for 
varying resolutions, the Xr profiles exacerbate these small discrep- 
ancies. At standard-res, Xr is an increasing function of radius, 
with 75% of all projections having values below 0.1 at 1 R e , and 
25% above 0.1. However, if we are taking into account the pro- 
jection which maximises Xr, the standard-res remnant should be 
classified as a fast-rotator. In the same context, both the high-res 
and very high-res are classified as slow-rotators. They have not 
the same profiles but have a similar overall behaviour: Xr first in- 
creases up to about 0.5 Re, and then decreases (up to 1.5 Re for 
the high-res and 1 Re for the very high-res). Outside 1.5 R e , Xr in- 
creases again outwards. Such a Ah profile is the clear s ignature of 
large-scale KDCs as m entioned in lEmsellem et al.li2007l (see also 
iMcDermid et all2006h . 
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Figure Al. Ellipticity and a±/a profiles (left and right panels resp.) of the four simulations as a function of R/ R e . From top to bottom : Simulations Dry-Fast, 
Wet-Fast, Dry-Slow and Wet-Slow. The three resolutions are shown with different colours : the standard-res in blue, the high-res in green and the very high-res 
in red. For each resolution, we plot five lines which correspond to the minimum and maximum at each radii (dashed lines), the mean value (thick solid lines) 
and the first and third quartiles (thin solid lines). The inter-quartiles space (which correspond to 50% of all projections) is filled with the colour associated with 
the resolution. 
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Figure A2. Left Panels The radial velocity and velocity dispersion profiles (in km s~ ) for the mean ellipticity projection along the global kinematic position 
angle (radius in kpc). Right Panels \r profiles as a function of R/R e . From top to bottom : Simulations Dry-Fast, Wet-Fast, Dry-Slow and Wet-Slow. The 
three resolutions are shown with different colours : the standard-res in blue, the high-res in green and the very high-res in red. 
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